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Abstract The rotation of stars has many interesting and important consequences 
for the photometric and chemical evolution of galaxies. Many of the predictions 
of models of stellar rotation are now compared with observations of surface abun- 
dances and velocities, with interferometric studies of fast rotating stars, with internal 
rotation profiles as they can be deduced by asteroseismology, to cite just a few obser- 
vational constraints. In this paper, we investigate how the outputs of models depend 
on the prescriptions used for the diffusion coefficients included in the shellular ro- 
tating models. After recalling the various prescriptions found in the literature, we 
discuss their impacts on the evolutionary tracks and lifetimes of the Main-Sequence 
(MS) phase, the changes of the surface composition and velocities during the MS 
phase, the distribution of the core helium lifetime in the blue and the red part of the 
HR diagram, the extensions of the blue loops, the evolution of the angular momen- 
tum of the core, and the synthesis of primary nitrogen in fast-rotating metal-poor 
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Impacts of various prescriptions 



massive stars. While some of these outputs depend only slightly on the prescrip- 
tions used (for instance, the evolution of the surface velocities), most of them show 
a significant dependence. The models which best fit the changes of the surface abun- 
dances are those computed with the vertical shear diffusion coefficient of Maeder 
(1997) and the horizontal shear diffusion coefficient by Zahn (1992). 



1 Rotation in stellar models 



In recent years, many effects of axial rotation on the structure and the evolution 
of ma ssive stars have been studied (see e.g. the recent review bv lMaeder & Mevnet 
120121) . Among the effects which are the most important are those linked to the trans- 
port of angular momentum and of chemical species in the interior of stars. These 
may strongly affect many outputs of stellar models such as the variation with the 
age of the surface abundances and velocities, the evolutionary tracks and lifetimes, 
the nature of the supernova events and of the stellar remnants, and the nature and 
the amount of new synthesized species. As a consequence, when results of rotating 
models are used in population synthesis models or in models for the chemical evolu- 
tion of galaxies, quite different results are obtained with respect to results obtained 
from non-rotating models. 

Most, if not all of the recent grids of rotating mo dels have been computed in the 
framework of the theor y proposed by IZahnl (119921) . with further improvements by 
Maeder & Zahnl dl998h . This was named the theory of shellular rotation, since it 



is based on the assumption that on an isobaric surface, the angular velocity, Q, is 
nearly constant, that means that any variations can be considered as a small pertur- 
bation. This nearly constant value of £2 on an isobaric surface is due to the fact that 
along those directions, there are neither stable temperature nor density gradients 
which counteract shear turbulence. This implies the existence of strong "horizon- 
tal" (i. e. along an isobaric surface) diffusion coefficient called Dh hereafter. In the 
following, when we speak about rotating models, we implicitly assume that we con- 
sider models with shellular rotatio n. The present models do not include the effects 
of the dynamo theory suggested by Sprui j (l2002lJ T l 

In the framework of the shellular theory of rotation, th e equation describing the 
transport of chemical species is a pure diffusive equation (IChabover & Zahnll 19921) 
written 
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where X, is the abundance in mass fraction of particles i, and D c k&m, the appropriate 
diffusion coefficient for chemical elements (see below). 



Numerical simulations bv lZahnet al] l l2007l) have studied MHD instabilities arising in the radia- 
tion zone of a differentially rotating star, in which a poloidal field of fossil origin is sheared into a 
toroidal field. Their simulations show no sign of dynamo action. 
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In a differentially rotating star, the evolution of the angular velocity Q has to be 
followed at each level r (for shellular rotation), so that a full description of Q(r, t) is 
available. The values of Q(r,t) influence the mixing of elements and in turn the evo- 
lution of Q (r, t) also depends on the mixing processes and on the distribution of the 
elements. The derivation of the equation for the transport of angular momentum is 
not straightfor ward. In the case of shellula r rotation, the equation in the Lagrangian 
form becomes dZahnll 1 992t lMaederll2009l) 

p|(^V,-^|:(P^^W) + ^(p^ g , 4 §) • (2) 

Here, Q is the average value of £2 on an isobar. U2 is the radial component of the 
meridional circulation velocity, and D ang , the appropriate diffusion coefficient for 
angular momentum. The second term on the right is a diffusion term, similar in its 
form to ((TJ, while the first term on the right is an advective term, i. e. modeling 
the transport by a velocity current. We notice that Eq. (U) does not contain such 
an advective term. It cou ld contain a term of that kind, however it can be shown 
dChaboyer&Zahnll992h that the combined effect of turbulence and circulation cur- 
rents is equivalent to a diffusion/or the element transport (see Eq.[T|l. 

In the equation for the transport of chemical species in radiative zones, the dif- 
fusion coefficient, D c h em , is made up of two terms. These are the vertical shear dif- 
fusion coefficient D s h ea r and the effective diffusion coefficient, Z) e ff , which account 
for the resultant effect of the strong horizontal shear diffusion, D n (i. e. the shear on 
an isobaric surface), and of the meridional currents. 

In the equation for the transport of angular momentum, the diffusion coefficient, 
Dang, is made up of only one term, the shear diffusion coefficient Z) s h ea r- 

For the coefficients £> s h e ar, we can find two different expressions in the literature. 

Z) shear from Maeder (1997, M97) 

Hp K (9n dlnX2\ 2 

Dshear " /energ ^[|V M + (V ad -V rad )] 1,32 ° dlnTj (3) 

where K = ^^f, and with / energ = 1, and <p = (0) ^ = 1. 
Dshear from Talon & Zahn (1997, TZ97) 



D s hear — fei 



Hp (K + D h ) (9% dlni2\ 



2 



- T , t\— ^2 (4) 

^[f V M(l + i)+(Vad-V rad )] 1 32 dlnr; 



with K, /energ, and <p as in (|T|>. 
For the coefficients Z) n , we can find three different expressions in the literature: 
D h from Zahn (1992, Z92) 



4 



Impacts of various prescriptions 



D h = -r\2V(r)-aU(r)\ (5) 

where a = \ d1 ^^ and Ch = 1 . 
D h from Maeder (2003, M03) 

D h =A r (r£2(r) V \2V-aU\) l,i (6) 
with a as in Eq. (5) and A = 0.002. 
D h from Mathis et al. (2004, MZ04) 

Dh= (^) 1/2 ^) 1/2(r|2y - af/|)1/2 (?) 

with a as in Eq. (5) and J3 = 1.5 • 10~ 6 . 
All prescriptions use the same effective mixing coefficient for the chemical species: 

D eff =lMC (8) 
ett 30 D h w 

There are therefore 6 different combinations of the two shear diffusion coeffi- 
cients and of the three horizontal diffusion coefficients. The physics sustaining the 
different expressions for these various diffusion coefficients is described in details in 
the papers indicated above and we shall not recall them in the present work. We just 
summarize below a few facts which are useful to keep in mind in order to understand 
their different impacts in stellar models. 

• Since the angular momentum is transported mainly by the meridional currents, 
one can expect that changing the expressions for the diffusion coefficients will 
have only a weak impact on the angular momentum distribution in stars. We 
shall see that this is well verified by the numerical models. 

• The diffusion coefficient D e ff is the key quantity determining the efficiency of 
mixing in regions where there is a strong /x-gradient, for instance at the border 
of the H-convective core. 

• The larger Dh is, the smaller will be D e ff, and thus less mixing will occur in 
regions of strong jj, -gradients. The expressions of M03 and MZ04 for Dh are 
larger than the expression given by Z92. 

• The diffusion coefficient D s h ea r is the key quantity determining the efficiency 
of mixing in regions with weak or no fx -gradients, typically in the radiative 
envelope of massive stars above the H-convective core. 

• The two expressions of D s h ea r are strictly equivalent in zones with no jii- 
gradients. 

• The ratio D s h ea r(M97)/D s hear(TZ97) <~ K/D^ in regions where ^V^ is signifi- 
cantly larger than the difference V a( j — V ra[ [. Since Dh is inferior to K, one has 

that D shea r(M97) > D she ar(TZ97). 
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In the present work, we study the implications of these different choices on the 
following model outputs: 

• Evolutionary tracks and lifetimes during the Main-Sequence (MS) phase. 

• Changes of the surface composition during the MS phase. 

• Changes of the surface velocities during the MS phase. 

• The distribution of the core helium lifetime in the blue and the red part of the 
HR diagram. 

• The extensions of the blue loops. 

• The evolution of the angular momentum of the core. 

• The synthesis of primary nitrogen in fast rotating metal poor massive stars. 

For that purpose we have computed models for different initial masses, metallicities, 
rotations with each one of the 6 possible combinations of values for (£> s hear, Ai) (see 
Table 1). For each mass and metallicity, the models are labeled by one digit and 
one letter: 1 is for models computed with the shear diffusion coefficient of Maeder 
(1997), and 2 for models computed with the shear diffusion coefficient of Talon 
& Zahn (1997), the letters A, B and C are respectively for the horizontal diffusion 
coefficient from Zahn (1992), Maeder (2003) and Mafhis et al. (2004). 

In Table 1, the first column gives the prescription used. The time-averaged equa- 
torial velocity during the MS phase is given in column 2, the MS lifetime is given 
in column 3, the difference between the surface helium abundance in mass frac- 
tion at the end of the MS phase and on the ZAMS is given in column 4. Column 
5 presents the N/H ratio obtained at the surface, at the end of the MS phase, and 
normalized to the initial N/H value. The core He-burning lifetime, and the analogs 
of columns 4 and 5 but at the end of the core He-burning phase are indicated in 
columns 6, 7 and 8 respectively. The duration of the core He-burning phase spent 
in the red (logoff < 3.68), in the blue (logoff > 3.87) and in the yellow part 
(3.68 < logoff < 3.87) of the HR diagram are given in columns 9, 10 and 11 re- 
spectively. The ratio of the time spent in the blue to that spent in the red is shown 
in column 12, the masses of the helium cores, of the carbon-oxygen cores and of 
the remnants are given in columns 13, 14 and 15. The mass of nitrogen produced 
divided by the mass of CNO elements initially present is given in column 16. 



3 Evolutionary tracks and lifetimes during the Main- Sequence 
phase 

Figure 1 presents the tracks in the Hertzsprung-Russell diagram (left panel) and the 
variation of the mass fraction of hydrogen as a function of the Lagrangian mass 
coordinate (right panel) for the 15 M at Z=0.002 and with i2; n i/i2 cr i t = 0.5. One 
can classify the models in three categories: 
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Table 1 Models computed in the present work 

iq 'h AY, to AY, ~ c ~i ~A |*J M He M C0 M ren , ^^Qj 

kms~' My My My My My M @ M . M e 

9 M: , Z=0.002, tta/Udi = 0-5 

1A 177 30.065 +0.0149 4.9 2.966 0.0157 4.9 3.0964 0.0970 0.0104 31.9 2.9 1.5 1.3 0.5 

IB -0.6% -11.8* +0.0051 3.8 +37.6% 0.0055 3.8 +37.5% +21.0% +27.9% 36.3 2.7 1.5 1.3 0.5 
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1 . The models which present very little differences with respect to the non-rotating 
model. These are the models 2B (not shown in Fig. 1) and 2C. 

2. The models which become overluminous with respect to the non-rotating model 
but present no significant extension towards lower effective temperatures. These 
are the models IA, IB (not shown) and 1C. 

3. Only one model, the model 2 A, becomes overluminous and reaches lower ef- 
fective temperatures at the end of the MS phase. 

This behavior reflects differences in the efficiency of mixing in different regions 
of the stars. The models which are of the less efficiently mixed stars (2B and 2C, 
due to large and thus small Z) e ff) show indeed very little difference in their 
non-rotating tracks (compare lines for the model 2C and the non-rotating model in 
Fig- [0- The model 2A presents a situation where there is an efficient mixing at the 
border of the convective core but where the shear is not so efficient in the radiative 
envelope. Thus the surface abundances are not yet modified at the stage represented 
in Fig.[TJ(right panel). The larger core increases the luminosity and also is responsi- 
ble for the extension towards lower temperature of the MS band in the HR diagram 
(an effect similar to an extension of the core produced by an overshoot). Finally 
the model which presents the greatest efficiency of mixing both at the border of the 
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Fig. 1 Evolutionary tracks in the Hertzsprung-Russell diagram (left panel) and the variation of the 
mass fraction of hydrogen as a function of the Lagrangian mass coordinate (right panel) for the 15 
M at Z=0.002 and with Q^/^cnt = 0.5. The dotted line corresponds to the non-rotating model, 
the other models are labeled as indicated in Table 1. The models used in the right panel have all a 
logTeff equal to 4.48. The central hydrogen mass fraction is equal to 0.1973 (non rotating model), 
0.1982 (1A), 0.2068 (1C), 0.2427 (2A), 0.2244 (2C). The models IB and 2B are not shown since 
they present many similarities with the models 1C and 2C. 



core and in the radiative zone is the model 1A. This mixing keeps this model in 
bluer positions in the HR diagram compared to model 2. 

One notes that in model 1 A, the /x-gradient (with respect to mass) is steeper than 
in the model 2A, while both models have the same expression for D e s (since they 
have the same expression for £>h). This is because in model 1A, which has £> s hear 
larger than in 2 A, hydrogen flows more efficiently inwards and helium outwards. 
This replenishes hydrogen at the border of the core. The net effect of these diffusions 
of hydrogen and helium is to make the star more luminous and bluer. The same 
occurs in the 1C model, although the effect is less marked because of the smaller 
value for D e ff. 

In Table 1, the MS lifetimes are indicated in column 3. The first row for each 
model gives the value of the MS lifetime in million years for the model 1A. The 
rows for the other models show the differences in percentage with respect to the 
value obtained for the model 1A. We see that the impact on the MS lifetimes re- 
mains modest (at most 13.5% for the 15 M© models considered here) with respect 
to the precision with which an age estimated can be made through the fitting of an 
isochrone in this mass domain. On the other hand the scatter is not negligible with 
respect to the amplitude of the effect of the increase of the MS lifetime due to ro- 
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tation. Indeed, the increase of the MS lifetime with respect to the non-rotating case 
amounts to 18.5% for the 15 M model. 

All the models computed here account for the same overshoot, but Fig.[T]shows 
that the track will present quite different extensions due to the various prescriptions. 
One sees also that the models 1A and 2 A (those computed with the smaller 
and thus greater D e g) are the only ones presenting an extension of the core with 
respect to the non-rotating model. Therefore only the use of these two prescriptions 
can attribute part of the extension of the convective core to an effect of rotational 
mixing. Let us note that the recent de t ermin ation of the extension of the mixed core 
in fast rotating stars by iNeiner et al.l (120121) seem to support the view that rotation 
enlarges the convective core. This would support prescriptions 1A or 2A. In that 
case, one should use slowly rotating stars in order to constrain the extension of the 
core due to the proce ss of convective penetration alone as is done, for instance, in 



Ekstrom et all [20 12). 



4 Changes of the surface composition during the MS phase 

The changes of the surface composition during the MS phases can be seen in Table 1 
and in Fig.|2]for the 15 M e at Z=0.002 and Q m i/Q cr ^ t = 0.5 cases. Column 3 and 
4 of Table 1 give, at the end of the MS phase, the excesses of helium at the surface 
of the star (in mass fraction and with respect to the initial value) and the ratio of 
nitrogen to hydrogen normalized to the initial value, respectively. 

Independent of the prescriptions used, o ne notes, as was already obtained in pre- 



vious works (see iMaeder & Mevnetll2001l) that the surface enrichment in nitrogen 



increases in increasing initial stellar masses and that it also increases for decreasing 
metallicity, in both cases keeping the initial rotation the same and comparing stars 
at similar evolutionary stages. 

One sees also that the surface enrichments are higher when £> s hear is higher (com- 
pare models of series 1 with models of series 2). This is quite logical since £> s hear 
is the parameter which governs the transport of chemical elements in the region 
extending from the vicinity of the core up to the surface. 

The greatest surface enrichments are always obtained for the model 1A, the 
smallest for the model 2A (except for the very metal poor 40 Mq model, but the 
difference between models 2A and 2B is quite small). For a given model, the vari- 
ation of the prescriptions used produces a scatter of the N/H value obtained at the 
end of the MS phase of around a factor of 2. For the excesses in helium, the factors 
are greater, however, except for the 40 M Q at Z=0. 00001 case, the enhancements are 
very modest and well below those which could be estimated from observed spectra. 
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Fig. 2 Upper panel Evolution of the N/H ratio at the surface of stellar models computed with 
various diffusive coefficients as a function of the equatorial velocity. Lower panel Evolution of the 
effective temperature during the core He-burning phases for different prescriptions of the diffusive 
coefficients. 
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5 Changes of the surface velocity during the MS phase 

Column 2 in Table 1 compares the time-averaged equatorial velocity during the MS 
phase. The first row of the table indicates the equatorial surface velocity for model 
1A in kms -1 . The other lines indicate the difference in percentage with respect to 
model 1A. One sees that the differences are always less than 5.3%, which is small. 
This is due to the fact that changing the prescriptions for the diffusion coefficients 
has only a very weak impact on the velocity of the meridional currents. Thus the 
transport of the angular moment which is mainly driven by these currents is almost 
unaffected by changes of the diffusion coefficients. 

This can also be seen in Fig. 2 (left panel), which shows the variation of the 
nitrogen enhancement with respect to the surface equatorial velocity (tracks on this 
diagram go from right to left as time proceeds). One sees that during the MS phase, 
all of the models span the same interval in equatorial velocities. What changes is the 
surface nitrogen enrichment obtained at the end of the MS phase. In the plane N/H 
versus u eq , model 1 A will produce a steeper relation than model 2 A. 

If one considers mean values of the ratio N/H for B dwarfs in the Galaxy and 
in the SMC one obtains enhancement fact ors between 1 .6 and 2.5 for the Galaxy 
and 3.2 and 6.3 for the SMC (see Table 2 in Ma eder & Mevnetl2012 , and references 



therein). These mean values should be obtained for stars with an average rotation 
during the MS phase of around 200 km s _1 . 

At solar metallicity, the 15 M Q model 1A well reproduces the observed enrich- 



ments ( Ekstrom et al . 2012). This is not a surprise, since the value of the parameter 
/enera in D s h ear , chosen equal to one, has been selected in order to fit these observed 
enrichments. With the same value of / energ , the 15 M e at Z=0.002 with an average 
rotation of 200 km s and prescriptions lApredicts an enhancement factor between 
3 and 5.3 in the last third of its MS lifetime, which is in the range of the observed 
values. In that case, no calibration has been made and the good fit supports this kind 
of model. The other prescriptions gives too low surface enrichments keeping / en erg 
equal to 1 . 

We note that a value equal to 1 for / en erg implies that we really account for the 
physics involved into the expression for D s h ear , which would not be the case if one 
would have to multiply the expression by a constant much greater or much smaller 
than one! 



6 The distribution of the core helium lifetime in the blue and the 
red part of the HR diagram 

The observed number of blue to red supergiants in clusters at different metallicities 
is an important feature that stellar models should be able to reproduce. This is im- 
portant for many reasons, for instance to predict the correct photometric evolution 
of young starburst regions or the nature of the progenitors of the core collapse su- 
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pernovae. It happens that observation shows that the blue to red ratio in clusters with 
masses at the turn off between 9 and 30 M Q increases when the metallicity increases, 
while standard stellar models pr edict that the blue to red s upergiant ratio decreases 
when the metallicity increases dMevlan & Maederi 1 19821; lLanger & Maederi 1 19951: 
lEggenberger et al.l l2002; Me vnet et al.ll201 lbl) . The blue to red s upergiant ratio ha s 
also been discussed in the contex t of field stellar populations in lHartwickl(ll970l) : 
iDohm-Palmer & Skillmanl d2002h . 

At the moment, there is no explanation for this general trend. On the other hands, 
many works could reproduce the blue to red supergia nt ratios observed at o ne given 
metal licity by changing the mass loss rates (see e.g. | Salasnich et al.lll999l) or mix- 
ing (lLanger & Maederi 1 19951 iMaeder & Mevnetll200ll) . In this paper we shall not 
discuss all the aspects of this question but focus on the importance of mixing. 

Looking at the right panel of Fig. [2] which shows the evolution of the effective 
temperature during the core Helium burning phase, we see that only one set of dif- 
fusion coefficients (the one using the D s h ea i- of Talon & Zahn 1997, and the of 
Zahn 1992) makes the 15 M at Z=0.002 evolve rapidly to the re d part of the HR di- 
agram after the MS phase. This is the prescription that we used in lMaeder & Mevnet 
(12.0011) . where we suggested that rotational mixing could help a lot in reproducing 
the observed blue to red supergiant ratio in the Small Magellanic Cloud. In view of 
the present results we see that, while this conclusion might always be correct, it is 
however quite dependent on the prescriptions used for the diffusion coefficients. 

It is interesting to identify from the numerical experiments performed in this 
work, the conditions which favor a rapid redward evolution at low metallicity. It 
does appear that two conditions have to be satisfied: 1) The mixing at the border of 
the convective cores (both during the H- and He-burning phases) have to be suffi- 
ciently efficient. Indeed, we see that any very strong values for Dh, which prevent 
any strong mixing in regions with strong molecular weight gradient (/x-gradient), 
also prevents the star evolving to the red phase. 2) The mixing in the zones where 
the /^-gradients are weak, namely the outer part of the radiative envelope, should not 
be too strong, because any strong mixing there would make the star be more "ho- 
mogeneous" and thus maintain a bluer position in the HR diagram. As very often 
noticed in the literature, we see that this red to blue evolution is a feature which is 
very sensitive to many physical ingredients of the models. The fact that it depends 
on subtle changes of the efficiency of mixing in different regions of the star is just 
one illustration of this. 

Taken at face value, the set of diffusion coefficients D s h ea r from Talon & Zahn 
(1997), and the Dj, of Zahn (1992) appear as the most favored to explain the blue 
to red supergiant ratio at low metallicity. However, other parameters — for instance 
the changes of the surface abundances expected for an averaged rotational velocity 
— are better fitted with the prescription 1A (keeping / en erg = 1)- Moreover, since 
mass loss, both during the MS phase and at the red supergiant phase, pl ays a key role 



in sha ping the blue to red supergiant ratio (see e.g. the discussion in IMevnet et al 



201 la) , it may be premature to use the observed variation of the blue to red su- 
pergiant ratio to constrain the prescription to be used. Probably in order to make 
progress in this area of research two important points have first to be settled: 1) to 
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distinguish, using observations of the surface abundances and/or of the vibrational 
properties, those blue supergiants which are direct succes sors of MS stars from those 
which are on a blue loop after a red supergiant stage. The ability to distinguish, at 
different metallicities, between those blue supergiants coming from the MS phase 
from those coming from the red supergiant phase would improve considerably our 
understanding of the blue supergiant formation process; 2) to obtain more reliable 
mass loss rates especially during the red and blue supergiant phases. 



7 The extensions of the blue loops 

Another feature which is sensitive to the form of the diffusion coefficients is the 
extension of the blue loops for stars with masses between about 3 to 12 M Q . This is 
important in order to predict the populations of Cepheids, and also for for the blue 
to red supergiant ratio discussed above, since the presence of a blue loop (when 
compared with the same model without blue loop) will reduce the lifetime of red 
supergiant but increase that of the blue supergiant. 

Looking at Fig.[5J we can see the following features: at Z = 0.002, the duration 
of the blue loops increases more and more passing from models 2A, to 2B, and then 
to 2C. When the D s h ea r is changed (model 1), the "loop" (if we can still speak of a 
loop in this case) even begins in the blue part of the HR diagram. 

At Z=0.014, the situation is quite different, first the loops in all models are sig- 
nificantly reduced, which is a well-known effect when the metallicity increases. 
Second, models 2B and 2C do not show any loops. In case such prescriptions would 
be adopted, then only slow rotators can show a blue loop and thus explain the exis- 
tence of Cepheids. If such prescriptions were to be adopted, then only slow rotators 
could show blue loops and thus explain the existence of Cepheids. 

We can note also that in the first set of models (labels beginning with one, i. e. 
D s hear from Maeder 1997), the impact of changing Dh on blue loops is quite modest 
in both metallicities. In the second set of models, changing Dh has a strong effect. 

At low metallicity, models 2A for 9 M e would be the more helpful to recon- 
cile the theoretical predictions with the observations as was suggested from previ- 
ous section which focused on 15 M stellar models. At solar metallicity, whatever 
model is considered, the blue to red supergiant ratios do appear too low (at most 1 
while it is observed at around 3). At this metallicity, the problem may be at least 
partially cured by enhancing the mass loss rate during the red supergiant phase. 



8 The angular momentum of the core 

During the evolution of a star, the core loses some angular momentum, mainly due 
to the effect of meridional currents. It happens that, in shellular rotating models 
without interior magnetic field, these losses are not sufficient to explain the rela- 
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Fig. 3 Upper panel Evolution of the effective temperature during the core He-burning phases for 
different prescriptions of the diffusive coefficients in a rotating 9 M Q stellar model at Z=0.002. 
Lower panel Same as the panel on the left at Z=0.014. 
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tively long observed rotation periods of young pulsars (IHeger et al.ll2.Q04l) . This is 
the reason why some authors have considered much stronger coupling between the 
core and the envelope by introducing a magnetic field wh ich forces solid b ody ro- 
tation or near solid body rotation during the MS phase (IHeger et al. 2005 ). Here, 
we have not accounted for such a strong coupling. A question however that we may 
ask is the extent to which this loss of angular momentum depends on the prescrip- 
tion used. Are there any of these prescriptions which would significantly change the 
angular momentum contained in the core at the end of its evolution? 

A priori, one would expect that the loss of angular momentum by the core due to 
the transport processes should be only slightly dependent on the various prescrip- 
tions because, as already stressed, during most of the stellar lifetime, angular mo- 
mentum is transported by the meridional currents whose velocities are only weakly 
dependent on the choice of Dh- Let us, however, check this point in the numerical 
models that we have. Since we have not pursued the computation beyond the end 
of the core He-burning phase, we compare here the angular momentum of the core 
obtained at the end of the core helium burning phase. The masses of the remnant of 
the different 15 M© Z=0.002 models are between 1.6 and 1.8 Mq. The angular mo- 
mentum which would be locked into the 1 .6 M© remnant supposing that no change 
does occur in the advanced phases of the evolution, would be between 0.93 and 1.14 
10 50 g cm 2 s _1 depending on the prescription used. The analog values in the case 
that the remnant is 1.8 M© would be between 1.19 and 1.44 1 50 g cm 2 s _1 . So we 
see that these quantities present a scatter around their mean values of at most 20%. 

Is such a scatter important? As it concerns the missing angular momentum loss 
of the core, the answer is clearly no. To illustrate this, let us derive the following 
numerical estimate: if we lock an angular momentum content of 1 10 50 g cm 2 s _1 in 
a neutron star, it would show a rotation period of about 0.1 ms, smaller by a factor 
between 4 and 7 than the critical periods for n eutron stars which are between 0.44 
and 0.65 ms as given by iGeorgy et al.l (12012b . The period is also smaller by two 
to three orders of magnit ude than the observed periods of young pulsar s which are 
between 20 and 100 ms dMushmov & Pagdfl996l: iMarshall et alj|l998h . We can at 
least conclude that the missing transport mechanism cannot be due to a particular 
choice of the diffusion coefficients for D s h ea r an d Dh, since whatever choice is made, 
the angular momentum content of the core is more or less the same at the end of the 
core He-burning phase. 

The angular momentum losses of the core may be underestimated either during 
the H and He-burning phases of the star and/or in the advanced phases and/or at the 
time of the supernova explosion and/or during the early years of the evolution of 
the new b orn neutron star. It may be that magnetic braking may play a role in this 
context dMeynet et alJl201 lah . 
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Fig. 4 Variation of the angular velocity as a function of the radius inside the 40 M,: . stellar models 
when the mass fraction of helium at the center is equal to 0.45. The dotted lines show, in each 
model, the variation as a function of the radius of the rate of production of energy by H-burning, 
£h, in units of 10 5 ergs g" 1 s _1 . 



9 The synthesis of primary nitrogen in fast-rotating metal-poor 
massive stars 

Fast-rotating massive stars may be t he sources of primary nitro gen in the early 
phases of the evolution of galaxies (iChiappini et al.ll200a 120081) . Therefore it is 
important to assess the extent to which the primary nitrogen production depends on 
the prescriptions used. In the last column of Table 1, we have indicated the mass of 
nitrogen (in solar masses) present in the region outside the stellar remnant normal- 
ized by the mass of CNO elements that were initially present in the same region of 
the star. We call this quantity M( 14 N)/M(CNO). When nitrogen is produced by the 
transformation of the carbon and oxygen initially present in the star (secondary ni- 
trogen production channel), then the quantity shown in Table 1 can at most be equal 
to one. It would be one if all the carbon and oxygen initially present in that region 
were to have been transformed into nitrogen. Actually it is less than one because, in 
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the regions inside the H-burning shell, the nitrogen which has been produced by the 
CNO burning will then have been further transformed into 22 Ne. 

We see that M( 14 N)/M(CNO) is inferior to one except in the fast-rotating very 
metal-poor 40 M P) stellar models. This i l lustrat es the results already discussed by 
Mev net & Maede j d2002h : iMevnet et alj J2006) that rotational mixing, by bring- 
ing carbon and oxygen (freshly synthesized in the helium-burning core) into the 
H-burning shell, enhances significantly the quantities of nitrogen produced. The 
quantities produced are no longer limited by the initial metallicity of the star, since 
the carbon and oxygen transformed into nitrogen are synthesized by the star itself 
through helium transformation (primary nitrogen production channel). 

We see, however, that the enhancements of primary nitrogen production present 
great variations depending on the prescriptions used for D s h ear and Dh. Before we 
analyze these results in more detail, let us recall a few general facts: 

1 . Primary nitrogen production depends on the efficiency of transport mechanisms 
in the region between the He-burning core and the H-burning shell; 

2. The diffusion coefficient which dominates the transport in that region is D s hear- 
A careful reader may be puzzled by such a statement since we mentioned above 
that D s hear operates mainly in weak /x-gradient regions while the appropriate 
diffusion coefficient in strong fi -gradient regions would be D e ff- Does it mean 
that the /x-gradients are not so strong at the border of the He-burning core? The 
answer is yes. Indeed, in the core-shell intermediate region, most of the time, 
the gradients of jj. are not very strong since the connected regions are all helium 
rich, moreover the gradient of Q which enters into the expression for D s h ea r is 
important. 

3. Since D shear is the dominant diffusion coefficient, the gradient of Q becomes 
the key factor for primary nitrogen production, a small gradient producing less 
efficient shear mixing than a steep gradient. 

In Fig. |U we can see the variations of £2 as a function of the radius in the 40 
M stellar models when the mass fraction of helium in the core is equal to 0.45 
(that means models at the middle of the core He-buming phase). We see that models 
using the prescription by Maeder (1997) have a smooth gradient immediately above 
the convective core and show in outer zones are a succession of regions with very 
steep and flat gradients. We shall call this zone the cliff. The cliff corresponds to the 
transition between the envelope and the core. The H-burning shell is more or less at 
its base (see Fig. 4). 

The smooth gradient results, at least in part, from the activity of the shear trans- 
port having occurred during the core He burning phase. The flat Q regions in the 
cliff are produced by intermediate convective zones which are no longer present at 
the stage shown in the figure but have appeared in previous stages of the evolution 
of the star. In the present models, we assume that convective regions rotate as solid 
bodies, hence the flattening of Q in these zones. The flat regions present a very low 
D s hear coefficient since the Q -gradients are very small. The part of the Q profile 
which is important for primary nitrogen production is the portion with a smooth 
gradient just above the He-burning core. 
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In the model IB (higher Dj, and thus smaller D e g), we note that the smooth 
gradient zone is more compact, making the gradient of Q. in that region steeper. 
This results from the less efficient chemical mixing at the border of the H-burning 
core. As a result, the helium core will be smaller and the transition zone between 
the core and the shell in a deeper part of the star. Since model IB presents a steeper 
gradient of Q. in the smooth gradient zone, it produces more primary nitrogen than 
model 1A (see Table 1). 

The configurations presented in the bottom panels, resulting from prescriptions 
2A and 2B, show different characteristics with respect to models 1A and IB in 
upper panels: 1) just above the core the gradient of Q. is steeper; 2) the region with 
a succession of strong i2-gradients and flat portions no longer exists. These two 
features result from less efficient shear transport at the border of the burning cores. 
Again, as in model 1, we see that when high values of are used, stronger gradients 
are obtained. 

We can see from Table 1 that, actually, the primary nitrogen production is larger 
in the models IB, 2 A and 2B than in model 1 A as can be expected from the line of 
reasoning above. An interesting conclusion of this discussion is that primary nitro- 
gen production depends mainly of the gradient of £2 just above the helium burning 
core. 



10 Conclusions 

We have studied the impact of various prescriptions on important outputs of rotating 
massive star models. The main conclusions are the following: 

• The outputs of stellar models which show marginal dependence on the prescrip- 
tions used for D s hear and Dj, are the MS lifetimes, the evolution of the surface 
velocities, and the evolution of the angular momentum of the core. 

• The outputs of stellar models which show significant dependence on the pre- 
scriptions used for D s h ear and Dh are the shape of the evolutionary tracks, the 
surface enrichments predicted for a given initial rotation, the blue to red evo- 
lution, and the extensions of blue loops and the amount of primary nitrogen 
produced. 

• The general trends of the increase of the mixing efficiency with the increase of 
the initial mass of the star, of its initial rotation and with the decrease of the 
initial metallicity remain the same whatever the prescriptions used for D s h ea r 
and Dh. 

The hope would be, of course, to identify the most realistic diffusion coefficients 
on the basis of physical considerations. From an analytic point of view, the difficulty 
is mainly in how to treat turbulence. This is reflected by the fact that in each of these 
expressions, a free parameter is present (/energ in D s h ea r f° r instance). At the moment 
this free parameter is chosen so as to allow models to reproduce one well identified 
observed feature (for instance, the surface enrichments observed at the surface of 
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MS B-type stars at solar metallicity). Thus a good fit with the observed feature 
which has been used to calibrate the models is not a test of the models. The tests 
of the models are comparisons with other observed features such as comparisons 
with surface enrichments for other velocities, metallicities or initial mass ranges. 
We show above that prescription 1A can account for the observed enrichments in 
the SMC, while the prescription has been calibrated on solar metallicity stars. Thus 
this can be viewed as a support of this kind of model. 

We see also that one set of prescriptions cannot give satisfactory fit to all the 
observed features discussed here. This is expected since these observed features are 
not all governed only by the way mixing is treated. As indicated in the paper, the 
blue to red supergiant ratio, for instance, also depends a lot on the way in which 
the mass loss due to stellar winds is implemented in the models. Moreov er, other 
effects, not accounted for here, such as close binary evolution, or magnetic braking, 
may also contribute to some of the observed features. 

In order to progress on the theoretical side, the key points reside in the capacity to 
treat turbulence in a more rigorous way. The key to progress on the theoretical front 
would seem to be the ability to treat turbulence in a more rigorous way. Probably 
multi dimensional hydrodynam ical simulations can provide important hints on this 
topic dArnett & Meakinll201 lh . 

Another complementary approach is to use well-observed stars to constrain the 
models. Stellar models which can account for many observed features have a great 
chance to provide a realistic description of the structure of stars. It may be some- 
times for wrong reasons, in the sense that the physical process invoked may not be 
exactly the one operating in nature, but the effect of this physical process, whatever 
it is, has some chance to produce the structure as obtained in the model which best 
fits the observations. 

We are confident that improvements in the two directions of hydrodynamical 
simulations and comparisons of evolutionary models with observed stars will allow 
us to constrain the possibilities for the physics occurring in stellar interiors. 
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manuscript. 
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